Structure of nascent 5S RNPs at the crossroad between ribosome assembly and MDM2-p53 pathways.

The 5S ribonucleoprotein (RNP) is assembled from its three components (5S rRNA, Rpl5/uL18 and Rpl11/uL5) before being incorporated into the pre-60S subunit. However, when ribosome synthesis is disturbed, a free 5S RNP can enter the MDM2-p53 pathway to regulate cell cycle and apoptotic signaling. Here we reconstitute and determine the cryo-electron microscopy structure of the conserved hexameric 5S RNP with fungal or human factors. This reveals how the nascent 5S rRNA associates with the initial nuclear import complex Syo1-uL18-uL5 and, upon further recruitment of the nucleolar factors Rpf2 and Rrs1, develops into the 5S RNP precursor that can assemble into the pre-ribosome. In addition, we elucidate the structure of another 5S RNP intermediate, carrying the human ubiquitin ligase Mdm2, which unravels how this enzyme can be sequestered from its target substrate p53. Our data provide molecular insight into how the 5S RNP can mediate between ribosome biogenesis and cell proliferation.

Molecular Dynamics Simulations of Reduced and Oxidized TFIIIA Zinc Fingers Free and Interacting with 5S RNA.

Interactions of zinc finger (ZF) proteins with nucleic acids and proteins play an important role in DNA transcription and repair, biochemical recognition, and protein regulation. The release of Zn2+ through oxidation of cysteine thiolates is associated with disruption of gene expression and DNA repair, preventing tumor growth. Multi-microsecond molecular dynamics (MD) simulations were carried out to examine the effect of Cys oxidation on the ZF456 fragment of transcription factor III A (TFIIIA) and its complex with 5S RNA. In the absence of 5S RNA, the reduced ZF456 peptide undergoes conformational changes in the secondary structure due to the reorientation of the intact ZF domains. Upon oxidation, the individual ZF domains unfold to various degrees, yielding a globular ZF456 peptide with ZF4 and ZF6, responsible for base-specific hydrogen bonds with 5S RNA, losing their ßßα-folds. ZF5, on the other hand, participates in nonspecific interactions through its α-helix that conditionally unravels early in the simulation. In the presence of RNA, oxidation of the ZF456 peptide disrupts the key hydrogen bonding interactions between ZF5/ZF6 and 5S RNA. However, interactions with ZF4 are dependent on the protonation state of His119.

Structural-profiling of low molecular weight RNAs by nanopore trapping/translocation using Mycobacterium smegmatis porin A.

Folding of RNA can produce elaborate tertiary structures, corresponding to their diverse roles in the regulation of biological activities. Direct observation of RNA structures at high resolution in their native form however remains a challenge. The large vestibule and the narrow constriction of a Mycobacterium smegmatis porin A (MspA) suggests a sensing mode called nanopore trapping/translocation, which clearly distinguishes between microRNA, small interfering RNA (siRNA), transfer RNA (tRNA) and 5 S ribosomal RNA (rRNA). To further profit from the acquired event characteristics, a custom machine learning algorithm is developed. Events from measurements with a mixture of RNA analytes can be automatically classified, reporting a general accuracy of ~93.4%. tRNAs, which possess a unique tertiary structure, report a highly distinguishable sensing feature, different from all other RNA types tested in this study. With this strategy, tRNAs from different sources are measured and a high structural conservation across different species is observed in single molecule.

Broken symmetry between RNA enantiomers in a crystal lattice.

Explaining the origin of the homochirality of biological molecules requires a mechanism of disrupting the natural equilibrium between enantiomers and amplifying the initial imbalance to significant levels. Authors of existing models have sought an explanation in the parity-breaking weak nuclear force, in some selectively acting external factor, or in random fluctuations that subsequently became amplified by an autocatalytic process. We have obtained crystals in which l- and d-enantiomers of short RNA duplexes assemble in an asymmetric manner. These enantiomers make different lattice contacts and have different exposures to water and metal ions present in the crystal. Apparently, asymmetry between enantiomers can arise upon their mutual interactions and then propagate via crystallization. Asymmetric racemic compounds are worth considering as possible factors in symmetry breaking and enantioenrichment that took place in the early biosphere.

Improving RNA Branching Predictions: Advances and Limitations.

Minimum free energy prediction of RNA secondary structures is based on the Nearest Neighbor Thermodynamics Model. While such predictions are typically good, the accuracy can vary widely even for short sequences, and the branching thermodynamics are an important factor in this variance. Recently, the simplest model for multiloop energetics-a linear function of the number of branches and unpaired nucleotides-was found to be the best. Subsequently, a parametric analysis demonstrated that per family accuracy can be improved by changing the weightings in this linear function. However, the extent of improvement was not known due to the ad hoc method used to find the new parameters. Here we develop a branch-and-bound algorithm that finds the set of optimal parameters with the highest average accuracy for a given set of sequences. Our analysis shows that the previous ad hoc parameters are nearly optimal for tRNA and 5S rRNA sequences on both training and testing sets. Moreover, cross-family improvement is possible but more difficult because competing parameter regions favor different families. The results also indicate that restricting the unpaired nucleotide penalty to small values is warranted. This reduction makes analyzing longer sequences using the present techniques more feasible.

Expansion segments in bacterial and archaeal 5S ribosomal RNAs.

The large ribosomal RNAs of eukaryotes frequently contain expansion sequences that add to the size of the rRNAs but do not affect their overall structural layout and are compatible with major ribosomal function as an mRNA translation machine. The expansion of prokaryotic ribosomal RNAs is much less explored. In order to obtain more insight into the structural variability of these conserved molecules, we herein report the results of a comprehensive search for the expansion sequences in prokaryotic 5S rRNAs. Overall, 89 expanded 5S rRNAs of 15 structural types were identified in 15 archaeal and 36 bacterial genomes. Expansion segments ranging in length from 13 to 109 residues were found to be distributed among 17 insertion sites. The strains harboring the expanded 5S rRNAs belong to the bacterial orders Clostridiales, Halanaerobiales, Thermoanaerobacterales, and Alteromonadales as well as the archael order Halobacterales When several copies of a 5S rRNA gene are present in a genome, the expanded versions may coexist with normal 5S rRNA genes. The insertion sequences are typically capable of forming extended helices, which do not seemingly interfere with folding of the conserved core. The expanded 5S rRNAs have largely been overlooked in 5S rRNA databases.

Evolutionary dynamics of 18S and 5S rDNA in autotriploid Carassius auratus.

The Carassius auratus (crucian carp) complex of the Dongting water system exhibits coexistence of diploid and triploid forms. As reported, triploid C. auratus is autotriploid origin. Ribosomal DNA (rDNA) with evolutionary conservation is widely used to study polyploidization. Here, we investigated genomic and transcribed rDNA sequences (18S and 5S) in diploid (2nCC, 2n = 100) and triploid (3nCC, 3n = 150) C. auratus. The results showed that the genetic traits and expression of 18S and 5S rDNA from 2nCC individuals were identified in 3nCC individuals. Moreover, pseudogenization of rDNA (18S and 5S) sequences were also observed in both 2nCC and 3nCC individuals, but expression of these variants was not detected. Based on the transcribed rDNA consensus sequence between 2nCC and 3nCC individuals, the functional secondary structures of 18S rRNA (expansion segments, ES6S) and 5S rRNA were predicted. These data demonstrated that complex evolutionary dynamics existed in the rDNA family of C. auratus. The evolutionary conservation of rDNA revealed that autotriploidization could not induce the divergence in Carassius taxa of the Dongting water system. These observations will expand our knowledge of the evolutionary dynamics of the rDNA family in vertebrates.

Mirror-Image 5S Ribonucleoprotein Complexes.

After realizing mirror-image genetic replication, transcription, and reverse transcription, the biggest challenge in establishing a mirror-image version of the central dogma is to build a mirror-image ribosome-based translation machine. Here, we chemically synthesized the natural and mirror-image versions of three ribosomal proteins (L5, L18, and L25) in the large subunit of the Escherichia coli ribosome with post-translational modifications. We show that the synthetic mirror-image proteins can fold in vitro despite limited efficiency and assemble with enzymatically transcribed mirror-image 5S ribosomal RNA into ribonucleoprotein complexes. In addition, the RNA-protein interactions are chiral-specific in that the mirror-image ribosomal proteins do not bind with natural 5S ribosomal RNA and vice versa. The synthesis and assembly of mirror-image 5S ribonucleoprotein complexes are important steps towards building a functional mirror-image ribosome.

RNA secondary structure prediction with pseudoknots: Contribution of algorithm versus energy model.

MOTIVATION: RNA is a biopolymer with various applications inside the cell and in biotechnology. Structure of an RNA molecule mainly determines its function and is essential to guide nanostructure design. Since experimental structure determination is time-consuming and expensive, accurate computational prediction of RNA structure is of great importance. Prediction of RNA secondary structure is relatively simpler than its tertiary structure and provides information about its tertiary structure, therefore, RNA secondary structure prediction has received attention in the past decades. Numerous methods with different folding approaches have been developed for RNA secondary structure prediction. While methods for prediction of RNA pseudoknot-free structure (structures with no crossing base pairs) have greatly improved in terms of their accuracy, methods for prediction of RNA pseudoknotted secondary structure (structures with crossing base pairs) still have room for improvement. A long-standing question for improving the prediction accuracy of RNA pseudoknotted secondary structure is whether to focus on the prediction algorithm or the underlying energy model, as there is a trade-off on computational cost of the prediction algorithm versus the generality of the method. RESULTS: The aim of this work is to argue when comparing different methods for RNA pseudoknotted structure prediction, the combination of algorithm and energy model should be considered and a method should not be considered superior or inferior to others if they do not use the same scoring model. We demonstrate that while the folding approach is important in structure prediction, it is not the only important factor in prediction accuracy of a given method as the underlying energy model is also as of great value. Therefore we encourage researchers to pay particular attention in comparing methods with different energy models.

Oxidative damage of 18S and 5S ribosomal RNA in digestive gland of mussels exposed to trace metals.

Numerous studies have shown the ability of trace metals to accumulate in marine organisms and cause oxidative stress that leads to perturbations in many important intracellular processes, including protein synthesis. This study is mainly focused on the exploration of structural changes, like base modifications, scissions, and conformational changes, caused in 18S and 5S ribosomal RNA (rRNA) isolated from the mussel Mytilus galloprovincialis exposed to 40µg/L Cu, 30µg/L Hg, or 100µg/L Cd, for 5 or 15days. 18S rRNA and 5S rRNA are components of the small and large ribosomal subunit, respectively, found in complex with ribosomal proteins, translation factors and other auxiliary components (metal ions, toxins etc). 18S rRNA plays crucial roles in all stages of protein synthesis, while 5S rRNA serves as a master signal transducer between several functional regions of 28S rRNA. Therefore, structural changes in these ribosomal constituents could affect the basic functions of ribosomes and hence the normal metabolism of cells. Especially, 18S rRNA along with ribosomal proteins forms the decoding centre that ensures the correct codon-anticodon pairing. As exemplified by ELISA, primer extension analysis and DMS footprinting analysis, each metal caused oxidative damage to rRNA, depending on the nature of metal ion and the duration of exposure. Interestingly, exposure of mussels to Cu or Hg caused structural alterations in 5S rRNA, localized in paired regions and within loops A, B, C, and E, leading to a continuous progressive loss of the 5S RNA structural integrity. In contrast, structural impairments of 5S rRNA in mussels exposed to Cd were accumulating for the initial 5days, and then progressively decreased to almost the normal level by day 15, probably due to the parallel elevation of metallothionein content that depletes the pools of free Cd. Regions of interest in 18S rRNA, such as the decoding centre, sites implicated in the binding of tRNAs (A- and P-sites) or translation factors, and areas related to translation fidelity, were found to undergo significant metal-induced conformational alterations, leading either to loosening of their structure or to more compact folding. These modifications were associated with parallel alterations in the translation process at multiple levels, a fact suggesting that structural perturbations in ribosomes, caused by metals, pose significant hurdles in translational efficiency and fidelity.

Advanced multi-loop algorithms for RNA secondary structure prediction reveal that the simplest model is best.

Algorithmic prediction of RNA secondary structure has been an area of active inquiry since the 1970s. Despite many innovations since then, our best techniques are not yet perfect. The workhorses of the RNA secondary structure prediction engine are recursions first described by Zuker and Stiegler in 1981. These have well understood caveats; a notable flaw is the ad-hoc treatment of multi-loops, also called helical-junctions, that persists today. While several advanced models for multi-loops have been proposed, it seems to have been assumed that incorporating them into the recursions would lead to intractability, and so no algorithms for these models exist. Some of these models include the classical model based on Jacobson-Stockmayer polymer theory, and another by Aalberts and Nadagopal that incorporates two-length-scale polymer physics. We have realized practical, tractable algorithms for each of these models. However, after implementing these algorithms, we found that no advanced model was better than the original, ad-hoc model used for multi-loops. While this is unexpected, it supports the praxis of the current model.

How Mg2+ ion and water network affect the stability and structure of non-Watson-Crick base pairs in E. coli loop E of 5S rRNA: a molecular dynamics and reference interaction site model (RISM) study.

The non-Watson-Crick (non-WC) base pairs of Escherichia coli loop E of 5S rRNA are stabilized by Mg2+ ions through water-mediated interaction. It is important to know the synergic role of Mg2+ and the water network surrounding Mg2+ in stabilizing the non-WC base pairs of RNA. For this purpose, free energy change of the system is calculated using molecular dynamics (MD) simulation as Mg2+ is pulled from RNA, which causes disturbance of the water network. It was found that Mg2+ remains hexahydrated unless it is close to or far from RNA. In the pentahydrated form, Mg2+ interacts directly with RNA. Water network has been identified by two complimentary methods; MD followed by a density-based clustering algorithm and three-dimensional-reference interaction site model. These two methods gave similar results. Identification of water network around Mg2+ and non-WC base pairs gives a clue to the strong effect of water network on the stability of this RNA. Based on sequence analysis of all Eubacteria 5s rRNA, we propose that hexahydrated Mg2+ is an integral part of this RNA and geometry of base pairs surrounding it adjust to accommodate the [Formula: see text]. Overall the findings from this work can help in understanding the basis of the complex structure and stability of RNA with non-WC base pairs.

Evidence of birth-and-death evolution of 5S rRNA gene in Channa species (Teleostei, Perciformes).

In higher eukaryotes, minor rDNA family codes for 5S rRNA that is arranged in tandem arrays and comprises of a highly conserved 120 bp long coding sequence with a variable non-transcribed spacer (NTS). Initially the 5S rDNA repeats are considered to be evolved by the process of concerted evolution. But some recent reports, including teleost fishes suggested that evolution of 5S rDNA repeat does not fit into the concerted evolution model and evolution of 5S rDNA family may be explained by a birth-and-death evolution model. In order to study the mode of evolution of 5S rDNA repeats in Perciformes fish species, nucleotide sequence and molecular organization of five species of genus Channa were analyzed in the present study. Molecular analyses revealed several variants of 5S rDNA repeats (four types of NTS) and networks created by a neighbor net algorithm for each type of sequences (I, II, III and IV) did not show a clear clustering in species specific manner. The stable secondary structure is predicted and upstream and downstream conserved regulatory elements were characterized. Sequence analyses also shown the presence of two putative pseudogenes in Channa marulius. Present study supported that 5S rDNA repeats in genus Channa were evolved under the process of birth-and-death.

Prediction of Secondary Structures Conserved in Multiple RNA Sequences.

RNA structure is conserved by evolution to a greater extent than sequence. Predicting the conserved structure for multiple homologous sequences can be much more accurate than predicting the structure for a single sequence. RNAstructure is a software package that includes the programs Dynalign, Multilign, TurboFold, and PARTS for predicting conserved RNA secondary structure. This chapter provides protocols for using these programs.

Involvement of human ribosomal proteins in nucleolar structure and p53-dependent nucleolar stress.

The nucleolus is a potent disease biomarker and a target in cancer therapy. Ribosome biogenesis is initiated in the nucleolus where most ribosomal (r-) proteins assemble onto precursor rRNAs. Here we systematically investigate how depletion of each of the 80 human r-proteins affects nucleolar structure, pre-rRNA processing, mature rRNA accumulation and p53 steady-state level. We developed an image-processing programme for qualitative and quantitative discrimination of normal from altered nucleolar morphology. Remarkably, we find that uL5 (formerly RPL11) and uL18 (RPL5) are the strongest contributors to nucleolar integrity. Together with the 5S rRNA, they form the late-assembling central protuberance on mature 60S subunits, and act as an Hdm2 trap and p53 stabilizer. Other major contributors to p53 homeostasis are also strictly late-assembling large subunit r-proteins essential to nucleolar structure. The identification of the r-proteins that specifically contribute to maintaining nucleolar structure and p53 steady-state level provides insights into fundamental aspects of cell and cancer biology.

Single TRAM domain RNA-binding proteins in Archaea: functional insight from Ctr3 from the Antarctic methanogen Methanococcoides burtonii.

TRAM domain proteins present in Archaea and Bacteria have a ß-barrel shape with anti-parallel ß-sheets that form a nucleic acid binding surface; a structure also present in cold shock proteins (Csps). Aside from protein structures, experimental data defining the function of TRAM domains is lacking. Here, we explore the possible functional properties of a single TRAM domain protein, Ctr3 (cold-responsive TRAM domain protein 3) from the Antarctic archaeon Methanococcoides burtonii that has increased abundance during low temperature growth. Ribonucleic acid (RNA) bound by Ctr3 in vitro was determined using RNA-seq. Ctr3-bound M. burtonii RNA with a preference for transfer (t)RNA and 5S ribosomal RNA, and a potential binding motif was identified. In tRNA, the motif represented the C loop; a region that is conserved in tRNA from all domains of life and appears to be solvent exposed, potentially providing access for Ctr3 to bind. Ctr3 and Csps are structurally similar and are both inferred to function in low temperature translation. The broad representation of single TRAM domain proteins within Archaea compared with their apparent absence in Bacteria, and scarcity of Csps in Archaea but prevalence in Bacteria, suggests they represent distinct evolutionary lineages of functionally equivalent RNA-binding proteins.

Simultaneous characterization of cellular RNA structure and function with in-cell SHAPE-Seq.

Many non-coding RNAs form structures that interact with cellular machinery to control gene expression. A central goal of molecular and synthetic biology is to uncover design principles linking RNA structure to function to understand and engineer this relationship. Here we report a simple, high-throughput method called in-cell SHAPE-Seq that combines in-cell probing of RNA structure with a measurement of gene expression to simultaneously characterize RNA structure and function in bacterial cells. We use in-cell SHAPE-Seq to study the structure-function relationship of two RNA mechanisms that regulate translation in Escherichia coli. We find that nucleotides that participate in RNA-RNA interactions are highly accessible when their binding partner is absent and that changes in RNA structure due to RNA-RNA interactions can be quantitatively correlated to changes in gene expression. We also characterize the cellular structures of three endogenously expressed non-coding RNAs: 5S rRNA, RNase P and the btuB riboswitch. Finally, a comparison between in-cell and in vitro folded RNA structures revealed remarkable similarities for synthetic RNAs, but significant differences for RNAs that participate in complex cellular interactions. Thus, in-cell SHAPE-Seq represents an easily approachable tool for biologists and engineers to uncover relationships between sequence, structure and function of RNAs in the cell.

5SRNAdb: an information resource for 5S ribosomal RNAs.

Ribosomal 5S RNA (5S rRNA) is the ubiquitous RNA component found in the large subunit of ribosomes in all known organisms. Due to its small size, abundance and evolutionary conservation 5S rRNA for many years now is used as a model molecule in studies on RNA structure, RNA-protein interactions and molecular phylogeny. 5SRNAdb (http://combio.pl/5srnadb/) is the first database that provides a high quality reference set of ribosomal 5S RNAs (5S rRNA) across three domains of life. Here, we give an overview of new developments in the database and associated web tools since 2002, including updates to database content, curation processes and user web interfaces.

Molecular organization of the 5S rDNA gene type II in elasmobranchs.

The 5S rDNA gene is a non-coding RNA that can be found in 2 copies (type I and type II) in bony and cartilaginous fish. Previous studies have pointed out that type II gene is a paralog derived from type I. We analyzed the molecular organization of 5S rDNA type II in elasmobranchs. Although the structure of the 5S rDNA is supposed to be highly conserved, our results show that the secondary structure in this group possesses some variability and is different than the consensus secondary structure. One of these differences in Selachii is an internal loop at nucleotides 7 and 112. These mutations observed in the transcribed region suggest an independent origin of the gene among Batoids and Selachii. All promoters were highly conserved with the exception of BoxA, possibly due to its affinity to polymerase III. This latter enzyme recognizes a dT4 sequence as stop signal, however in Rajiformes this signal was doubled in length to dT8. This could be an adaptation toward a higher efficiency in the termination process. Our results suggest that there is no TATA box in elasmobranchs in the NTS region. We also provide some evidence suggesting that the complexity of the microsatellites present in the NTS region play an important role in the 5S rRNA gene since it is significantly correlated with the length of the NTS.